CLAIMS 



1. A Voigt function approximation calculating method used for line-by-line 

calculations, said method comprising: 

(1) a step of dividing a domain of a Voigt function into a first range around the peak 
of the Voigt function and a skirt portion not contained in the first range, replacing the 
first range with a cubic function, calculating the values and derivatives of said cubic 
function and the Voigt function in the skirt portion for each of first predetermined 
intervals, and connecting said cubic function and said Voigt function at points of 
connection thereof using the values and derivatives of both functions; 

(2) a step of adding together the results of step (1) for a plurality of absorption lines; 

(3) a step of calculating the function values and derivatives for the results of step (2) 
by interpolation over intervals smaller than said first predetermined intervals; 

(4) a step of dividing said first range into a second range near the peak and a skirt 
portion not contained in the second range, replacing said second range of a "function 
representing the difference between the Voigt function and said cubic function" with a 
cubic function, and calculating values and derivatives of said cubic function and said 
"function representing the difference between the Voigt function and said cubic 
function" in the skirt portion for each of second predetermined intervals; 

(5) a step of connecting said cubic function and said "function representing the 
difference between the Voigt function and said cubic function" at points of connection 
thereof using the values and derivatives of both functions; 

(6) a step of adding the results of steps (4) and (5) to the results of step (3) for a 
plurality of absorption lines; 

(7) a step of calculating the function values and derivatives for the results of step (6) 

-15- 



by interpolation over intervals smaller than said second predetermined intervals; and 
(8) a step of adding the values of the "function representing the difference between 
the Voigt function and said cubic function" to the results of step (7) for a plurality of 
absorption lines in said second range. 

2. A method in accordance with claim 1, wherein (9) a step of calculating function 
values and derivative values by interpolation is repeated while narrowing the intervals until 
they become intervals of a minimum unit. 

3. A method in accordance with claim 1, wherein said steps (4) through (7) are 
repeated one or more times until the third predetermined interval is obtained. 

4. A method in accordance with claim 1, wherein said predetermined intervals are 
determined by using the following equation. 

The first predetermined interval for the widest sub-function is j kmax dv. Here, is a 
single-digit natural number, dv is the increment in wave number, and kmax is the largest 
natural number satisfying the relationship j kmax+2 pdv < Vmax. However, Vmax represents the 
maximum calculation range from the center of the absorption line, and p is a natural 
number for controlling the calculation precision 

5. A method in accordance with claim 3, wherein said predetermined intervals are 
determined using the following equation. 

The most specific third predetermined interval is j kmin dv. Here, ; is a single-digit 
natural number, dv is the increment of the wave number, and kmin is the maximum 
non-negative decimal fraction satisfying the relationship j kmin pdv < a (a being approximately 
y/4). However, y is an approximate value of the full-width at half-max of the absorption 
line, and p is a natural number for controlling the calculation precision. 
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6. A method in accordance with claim 1, wherein said second predetermined 

intervals are determined by using the following equation. 

For the sub-function with the (k - kmin + l)-th smallest width, the predetermined 
interval is j k dv. Here, is a single-digit natural number, dv is the increase in wave number 
and k is such that kmin <k< kmax. 



7. A method in accordance with any one of claims 1-3, wherein said interpolation is 

calculated with / set to 4, using the function values yo, yi and function derivative values yo', 
yi at xo, xi in the interpolation interval (xo, xi), with the below-given Equation (1) as a 
function value interpolation equation, Equation (2) as a function derivative value 
interpolation equation and e as a non-negative decimal fraction. 
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8. A method in accordance with any one of claims 1-3, wherein said interpolation is 

calculated with ; set to 5, using the function values yo, yi and function derivative values yo, 
yi' at xo, xi in the interpolation interval (xo, xi), with the below-given Equation (3) as a 
function value interpolation equation, Equation (4) as a function derivative value 
interpolation equation and e as a non-negative decimal fraction. 
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9. A method in accordance with any one of claims 1-8, for achieving high speeds, 

when assuming the Voigt function to be K(x, y) and the difference from the Voigt profile of 
the absorption line to be K(x, y) +/(*), replacing 

K(x,y) 
with: 

K(x,y)=AK(x,y)+Bf(x) 

and 

dK(x,y) 
dx 



with 



dK(x 9 y) = A dK{x 9 y) x B df{x) 

dx dx dx 



10. A method in accordance with any one of claims 1-8, for achieving high speeds, 

when assuming the Voigt function to be K(x, y) and the difference from the Voigt profile of 
the absorption line to be K(x, y)f(x), replacing 

K(x,y) 
with: 
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K(x,y) = K(x,y)f(x) 



and 



dK(x,y) 
dx 



with 



dK^y) = dK^y) () ( ^ 

dx dx J w v y} dx 



11. A method in accordance with any one of claims 1-10, for achieving high speeds in 

a sub-Lorentzian correction, by using 

K(x,y) = K(x,y)Aexp(-B\x\) 

and 

dK^y) = dK^y) ^ exp( _ ^ )+ ^ ^_ ^ ^ exp( _ ^ }] 



dx 



12. A method in accordance with any one of claims 1-10, for achieving high speeds in 

line-mixing correction, by replacing 

K(x,y) 
with: 

K(x,y) = AK(x,y)+BL(x,y) 

and 

dK(x,y) 
dx 



with 

dK{x,y) 
dx 



= -2 



{Ax + By)K{x,y)- {Ay - Bx)L{x,y)- 



B 



Here, L(x, y) is the imaginary component of the function w(z) (the real part is the Voigt 
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function), where complex number z = x + iy, defined by the following equation. 
(erfc(z) is a complex complementary error function). 

13. A Voigt function approximation calculating program used for line-by-line 

calculations, said program performing: 

(1) a step of dividing a domain of a Voigt function into a first range around the peak 
of the Voigt function and a skirt portion not contained in the first range, replacing the 
first range with a cubic function, calculating the values and derivatives of said cubic 
function and the Voigt function in the skirt portion for each of a first predetermined 
intervals, and connecting said cubic function and said Voigt function at points of 
connection thereof using the values and derivatives of both functions; 

(2) a step of adding together the results of step (1) for a plurality of absorption lines; 

(3) a step of calculating the function values and derivatives for the results of step (2) 
by interpolation over intervals smaller than said first predetermined intervals; 

(4) a step of dividing said first range into a second range near the peak and a skirt 
portion not contained in the second range, replacing said second range of a "function 
representing the difference between the Voigt function and said cubic function" with a 
cubic function, and calculating values and derivatives of said cubic function and said 
"function representing the difference between the Voigt function and said cubic 
function" in the skirt portion for each of second predetermined intervals; 

(5) a step of connecting said cubic function and said "function representing the 
difference between the Voigt function and said cubic function" at points of connection 
thereof using the values and derivatives of both functions; 

(6) a step of adding the results of steps (4) and (5) to the results of step (3) for a 
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plurality of absorption lines; 

(7) a step of calculating the function values and derivatives for the results of step (6) 
by interpolation over intervals smaller than said second predetermined intervals; and 

(8) a step of adding the values of the "function representing the difference between 
the Voigt function and said cubic function" to the results of step (7) for a plurality of 
absorption lines in said second range. 

14. A program in accordance with claim 13, wherein a step of calculating function 
values and derivative values by interpolation is repeated while narrowing the intervals until 
they become intervals of a minimum unit. 

15. A program in accordance with claim 13, wherein said steps (4) through (7) are 
repeated one or more times until the third predetermined interval is obtained. 
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